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Using lattice methods in non-canonical quantum statistics 
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We define a natural coarse-graining procedure which can be applied to any closed equilibrium quantum system 
described by a density matrix ensemble and we show how the coarse-graining leads to the Gaussian and canonical 
ensembles. After this motivation, we present two ways of evaluating the Gaussian expectation values with lattice 
simulations. The first one is computationally demanding but general, whereas the second employs only canonical 
expectation values but it is applicable only for systems which are almost thermodynamical. 



1. Introduction 

The usefulness of the canonical ensemble in sta- 
tistical mechanics is remarkable. The standard 
explanation of this success relies in taking the 
"thermodynamical limit" which corresponds to 
increasing the volume of the system to infinity 
while keeping all the relevant intensive quantities, 

1. e. densities, fixed and finite. From this point 
of view, the canonical ensemble should not have 
much utility for small systems consisting of only 
a few particles. This, however, does not seem 
to be the case, and in the following we propose 
a new approach which explains why and how the 
canonical ensemble can help also in the analysis of 
small equilibrium systems. We also explain how 
lattice simulations in general can be employed in 
this analysis. 

2. Ensembles from coarse-graining of the 
energy fluctuation spectrum 

The standard approach to quantum statistics 
Q uses a density matrix p, which is a non- 
negative, hermitian, trace class operator normal- 
ized to one and which gives the expectation value 
of an observable A by the formula (A) = Tr (Ap) . 
In some complete eigenbasis \ipi), the density ma- 
trix can thus be expanded as p — J2iPi\ipi)(ipi\, 
where the eigenvalues pi satisfy pi > and the 
normalization condition = 1- 

Suppose now that the system has a discrete 



energy spectrum, which in quantum mechanics 
is achieved for every potential that grows suffi- 
ciently fast at infinity. An equilibrium, i.e. time- 
independent, ensemble is then given by a density 
matrix which has time-independent eigenvalues 
and which satisfies [H, p] = 0. In this case, the 
eigenvectors ipi can be chosen so that they are 
also energy eigenvectors with an eigenvalue Ei. 

After these preliminaries, it is not hard to see 
that for any equilibrium ensemble which is given 
by a density matrix and which has energy as the 
only relevant conserved quantity, we can find a 
smooth mapping F so that p = F(H), i.e. that 
p = J2i F(Ei)\il)i){il)i\. We will call such a smooth 
mapping the fluctuation spectrum of the ensemble 
and we will use the term precanonical ensemble 
for those equilibrium ensembles which satisfy: 

1. The canonical partition function converges, 
Tre -0 < 00) f or a ii p > o. 

2. The energy fluctuations decay at least ex- 
ponentially at high energies: e l3E F(E) is a 
rapidly decreasing function ||] for all (3 < 

where /3+ > is a parameter. 

The following representation is then valid for 
any precanonical ensemble and for all < f3 < /3 + 
as long as A is a positive observable which satisfies 
Tr [Ae~@ H ] < oo in the same range of /3: 
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where F is the Laplace transform of F and the in- 
tegrand in the above equation is an analytic func- 
tion in the half-plane < Kew < f3+. This re- 
sult follows from Fourier-transform formulae for 
rapidly decreasing functions; the precise mathe- 
matical details can be found from ||]. 

The value of the integral representation ([!]) can 
be computed by saddle point methods and, for 
instance when A = 1, there is a unique positive 
saddle point (3 which dominates the value of the 
integral. For A = 1 it can be solved for small j3 
from the saddle point equation 

(H} c an = E + I3e 2 + 0(/3 2 ), (2) 

where, assuming the normalization JF = 1, 

E = JdxF(x)x, e 2 = JdxF(x)(x - E) 2 . (3) 

Since the precise form of the fluctuation spec- 
trum is difficult to measure, we need some way 
of parameterizing its large scale properties. As 
usual, these can be extracted from the orig- 
inal fluctuation spectrum by a coarse-graining 
transformation — here we used a convolution with 
a Gaussian distribution, 

/OO -I 
dyF{y)^= G ^ x '^ . 

Under this transformation, the Laplace-transform 
in (g) will change to F A (w) = ei A2w2 F(w). 

When A approaches infinity it is clear that the 
positive saddle point value (3 must go to zero. 
Since Fa is analytic near the origin, we can in 
this limit use the approximation Ew + ^e 2 w 2 
for InF^w) — the parameters are obtained from 
(||) by replacing F with Fa. Taking the in- 
verse Laplace-transform then shows that this cor- 
responds to using the Gaussian ansatz F\(x) = 
G e (E — x) for the fluctuation spectrum. This 
should not come as a surprise; the argumenta- 
tion is the same as used with the central limit 
theorem of probability theory. 

For large A the positive saddle point typically 
becomes dominant. On the other hand, the trace 
left in the positive saddle point approximation is 
simply the canonical trace, Tr (Ae~^ H ) : and often 
the canonical expectation value becomes a good 



approximation of the coarse-grained one. The 
precise condition for the use of the canonical en- 
semble can be given in terms of the canonical vari- 
ance a 2 = ((H — {H)) 2 ) c p u and the normalized 

canonical energy operator h = (H — (H))/a. The 
condition for using the canonical approximation 
for the partition function and the one for using 
the canonical expectation value for a positive ob- 
servable A are, respectively, 

2 

a = — r < 1 and a (Ah 2 ) I (A) < 1. 

Similarly, the Gaussian ensemble can be used if 
the left hand sides in the previous equations are 
not too large. A more complete explanation of 
these results can be found in || . 

The canonical approximation of the Gaussian 
expectation values has already been analysed in 
Q and we quote here only the results: The simple 
bounds, already referred to in the above, for the 
approximation of positive observables are 

(ih 2 ) ^ , (a)t: s ^ 

—a — — — < In — — < a 

(A) (A)f* 

and this approximation can be improved for awl 
by using the asymptotic series 

(A)f™ s = (A)^ D + ((l-h 2 )A)a + 0(a 2 ). (4) 

3. Gaussian ensemble on a lattice 

There are two different ways of using lattice 
simulations in the evaluation of Gaussian expec- 
tation values. In the direct approach, the lat- 
tice approximation is applied to the complex tem- 
perature trace in (|l|) which, after an exchange 
of the order of the integration and the contin- 
uum limit, leads to an integral kernel for the lat- 
tice simulations. Unfortunately the kernel is not 
a positive function and the results are in most 
cases obtained as a delicate cancellation of os- 
cillations of the kernel. This, however, is likely 
to be an unavoidable feature of any space-lattice 
simulation of microcanonical expectation values 
as the energy wavefunctions themselves are typi- 
cally highly oscillatory. 

This approach is most useful when canonical 
simulations at a complex temperature are pos- 
sible. Then the trace in the the integrand in (|l|) 
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can be evaluated in a number of points and the in- 
tegral computed by a discrete Fourier-transform. 
This would yield results for a range of values of 
the parameters E and e and, therefore, it would 
enable an inspection of a whole energy range at 
the same time. The main difficulty in this ap- 
proach is, of course, in the complex temperature 
lattice simulation with its oscillatory kernel func- 
tion. 

If an expectation value at only one value of 
the parameters E and e is needed, then a sec- 
ond alternative is also possible: perform first 
the Fourier-transformation of the canonical lat- 
tice kernel and do the lattice simulations with the 
resulting kernel function. In this case, the evalu- 
ation of the kernel function becomes an obstacle 
slowing down the simulation. 

The second approach to the Gaussian evalua- 
tion problem is to use the asymptotic series given 
in (Q) which requires the computation of canoni- 
cal expectation values only. The problem in this 
case is to find the correct lattice operators which 
would correspond to the different powers of h 
in the continuum limit. We will now show how 
these can be found in a simple quantum mechan- 
ical case and comment on some general features 
which should be relevant also for field theory lat- 
tice simulations — a more complete analysis of this 
kind of lattice system can be found from || . 

Consider, for simplicity, a non-relativistic par- 
ticle in a potential V{x). The Hamiltonian of this 
system is H = -^p 2 + V(x) and if the potential 
is bounded from below and increases sufficiently 
fast at infinity, the complex temperature trace has 
a rigorous lattice approximation given by 

Tre~ w " = lim U l x[L/{2ttw)] l / 2 C -^ Pi '- wVl , 

L^co J 
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A straightforward differentiation of this result 
then gives the operators which will measure in the 
continuum limit the expectation values of differ- 
ent powers of the Hamiltonian. We have given the 
first few of them in Table [l]. Two features of these 
results are worth pointing out: first, the kinetic 
energy is given by the operator \L — jjPl, which 



Table 1 

Lattice energy operators in terms of M 
P = \P L and Vp = PV L . 



L/2, 



"/3H" = ci ci 
"/3 2 # 2 " = c 2 + c 2 c 2 

«p3 H 3» =c 3 + 3ciC2 + C3 C3 



M -P f j + V f j 
M - 2P 
2{M-3P ) 



shows that Pl diverges as the lattice size in the 
continuum limit and thus needs to be "renormal- 
ized" . This reflects the well-known result that the 
continuum path-integral is concentrated on paths 
which are continuous, but non-differentiable. Sec- 
ondly, each power of the Hamiltonian needs a sep- 
arate renormalization term in the sense that using 
the powers of the operator giving the expectation 
value, ci, is not enough. This is exactly analo- 
gous to the situation of composite operators in 
field theory. 

4. Conclusions 

We have introduced the Gaussian ensemble as 
a means of refining the accuracy of the canon- 
ical ensemble and we have shown by a coarse- 
graining procedure why this would have applica- 
tions also for non-thermal equilibrium systems. 
The canonical, complex temperature, lattice sim- 
ulations offer one way of inspecting the behav- 
ior of the Gaussian expectation values. A second 
way, applicable for systems near the thermody- 
namical limit, uses correction terms which can be 
computed with the well-established methods of 
canonical lattice simulations. 
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